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ABSTRACT: The effect of rotation on the Rayleigh-Benard convection in a horizontal layer of ferrofluid is investigated 
by using Gale rk in weighted residual:, meiie.nl. Linear stability theory based upon normal made analysis and perturbation 
method is used to find expres \ions for Rayleigh number for free-free boundary layer of fluid. It is observed that the system is 
more stable in the rotating fluid than non-rotating fluid layer. 'Principle of exchange of stabilities' not valid and the 
oscillatory convection is possible only for certain conditions. ect of rotation and magnetic parameters on the 

stationary convection is investigated analytically. 
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Nomenclature 



a wave number 

B magnetic induction 

g acceleration due to gravity 

H magnetic field Intensity 

k thermal conductivity 

Ki pyomagnetic coefficient 

M magnetization 

Mi buoyancy magnetization 

M 3 magnetic parameter 

N magnetic thermal Rayleigh number 

n growth rate of disturbances 

p pressure (Pa) 

P r Prandtl number 

q fluid velocity 

R Rayleigh number 

t time 

T temperature 

Ta average temperature 

T A Taylor number 

u, v, w fluid velocity components 

(X, y, z) space co-ordinates 

Greek symbols 

a thermal expansion coefficient 

P uniform temperature gradient 

|i 0 magnetic permeability 

u viscosity 

p density of the fluid 

(p c ) heat capacity of fluid 

K thermal diffusivity 

O' perturbed magnetic potential 

co dimensional frequency 

% magnetic susceptibility 

Superscripts 

non dimensional variables 
perturbed quantity 

Subscripts 

0 lower boundary 

1 upper boundary 
H horizontal plane 
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I. INTRODUCTION 

Ferromagnetic fluid has wide ranges of applications in instrumentation, lubrication, printing, vacuum technology, 
vibration damping, metals recovery, acoustics and medicine, its commercial usage includes vacuum feed through for 
semiconductor manufacturing in liquid-cooled loudspeakers and computer disk drives etc. Owing the applications of the 
ferrofluid its study is important to the researchers. A detailed account on the stability of ferrofluid has been given by 
Rosensweig (1985) in his monograph. This monograph reviews several applications of heat transfer through ferrofluid. One 
such phenomenon is enhanced convective cooling having a temperature -dependent magnetic moment due to magnetization 
of the fluid. This magnetization, in general, is a function of the magnetic field, temperature, salinity and density of the fluid. 

In our analysis, we assume that the magnetization is aligned with the magnetic field. Convective instability of a 
ferromagnetic fluid for a fluid layer heated from below in the presence of uniform vertical magnetic field has been 
considered by Finlayson (1970). He explained the concept of thermo-mechanical interaction in ferromagnetic fluids. 

Thermo convective stability of ferromagnetic fluids without considering buoyancy effects has been investigated by 
Lalas and Carmi (1971). Linear and nonlinear convective instability of a ferromagnetic fluid for a fluid layer heated from 
below under various assumptions is studied by many authors Shliomis (2002), Blennerhassett et.al.(1991), Gupta and Gupta 
(1979), Stiles and Kagan (1990), Sunil et.al. (2005, 2006), Sunil, Mahajan (2008), Venkatasubramanian and Kaloni (1994), 
Zebib (1996), Mahajan (2010). Rotation also play important role in the thermal instability of fluid layer and has applications 
in rotating machineries such as nuclear reactors, petroleum industry bio mechanics etc. Owing to the various applications 
of ferrofluid an attempt has been made to investigate the thermal instability of a ferromagnetic fluid in the presence of 
rotation in a fluid layer heated from below using Galerkin weighted residuals method. 



II. MATHEMATICAL FORMULATION OF THE PROBLEM 

Consider an infinite, horizontal layer of an electrically non-conducting incompressible ferromagnetic fluid of 
thickness 'd', bounded by plane z = 0 and z = d. Fluid layer is acted upon by gravity force g (0, 0, -g) and a uniform 

magnetic field H = H^'k acts outside the fluid layer. The layer is heated from below such that a uniform temperature 



gradient (3| 
Fig.l 
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is to be maintained. The temperature T at z = 0 taken to be T 0 and T[ at z = d, (T 0 > T, as shown in 



: CO, O, -aO 
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Fig.l Geometrical configuration of the problem 

The governing equations under Boussinesq approximation for the above model (Finlayson (1970), Resenweig 
(1997), and Mahajan (2010) are: 

V.q = 0, (1) 

Po ^ = "Vp + p 0 g + |^V 2 q + n 0 (M.V)H + 2p 0 (q x Q) , (2) 
at 

( Po C 0 ) f ^ + (p o C 0 ) f q.VT = kV 2 T, (3) 
Maxwell's equations, in magnetostatic limit: 

V.B = 0,VxH = 0, B = n 0 (H + M). (4) 

The magnetization has the relationship 

M = |[M 0 +x(H-H 0 )-K 1 (T-T a )]. (5) 
The density equation of state is taken as 

p = p„[l-a(T-T a )]. (6) 

Here p, p 0 , q, t, p, \x, \x 0 , H, B, C 0 , T, M, K b and a are the fluid density, reference density, velocity, time, pressure, dynamic 
viscosity (constant), magnetic permeability, magnetic field, magnetic induction, specific heat at constant pressure, 
temperature, magnetization, thermal conductivity and thermal expansion coefficient, T a is the average temperature given by, 
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T a = j H = |H| , M = |M| and M 0 = M(H 0 ,T a ). The magnetic susceptibili 

and K 

UHj Ho , Ta 



ity and pyomagnetic coefficient 



„ ( dM\ 

are defined by % = and K, = respectively. 

\9T y Ho>Ta 

Since the fluid under consideration is confined between two horizontal planes z = 0 and z = d, on these two planes 
certain boundary conditions must be satisfied. We take case of free-free surface and assume the temperature and volumetric 
fraction of the nanoparticles are constant and the boundary conditions [Chandrasekhar (1961)] are 

w = 0, T = T 0 at z = 0, and w = 0, T = T, at z = d. (7) 
II. 1 Basic Solutions 

The basic state is assumed to be a quiescent state and is given by 

q(u,v,w) = q b (u,v,w) = 0 p = p b (z) T = T b (z) = -fiz + T a J K,(T b -T a )' 



1 + X 



K 2 (T b -T a 
1 + X 



, H o +M o =H o ext. 





, fM 0 X 











II.2 The Perturbation Equations 

We shall analyze the stability of the basic state by introducing the following perturbations: 

q = q b +q\ p = p b (z)+5p,T = T b (z) + 6 , H = H b (z)+H' M = M b (z)+ M' . (9) 
where q'(u,v,w), 5p, 0, H'(H'i,H' 2 ,H'3) and M'(M'i,M'2,M' 3 ) are perturbations in velocity, pressure, temperature, magnetic 
field and magnetization. These perturbations are assumed to be small and then the linearized perturbation equations are 
V.q'=0, (10) 

p o ^ = -V8p + H V 2 q' + Po ag6k - f(l + X )^k - K,ek jl + 2 Po (q x Q), (1 1) 

a i+x v Sz ') 

— = kV 2 9 + Pw, (12) 
dt 

^ = K,* (13) 

dz 2 dz 

Where H' = VcpJ and (p' is the perturbed magnetic potential and K = is thermal diffusivity of the fluid. 

(Po c o) f 

And boundary conditions ai'e 

w = 0, T = T 0 ,Dcp = 0 at z = Oand w=0, T = T p D(p = 0 at z = d. (14) 
We inttoduce non-dimensional variables as 

There after dropping the dashes ( " ) for simplicity. 

Equations (10) - (14), in non dimensional form can be written as 

V.q = 0, (15) 

— — = -V5p + V 2 q + R(l + M 1 )9k-RM 1 ^k + JT7(ve x -uej, (16) 
Pr dt dz 

— = V^6 + w • (17) 
at 

M 3 V 2 <j, 1 -(M 3 -l)^% = ^. (18) 

dz dz 
Where non-dimensional parameters are: 
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p _ J£^ is Prandtl number; R _ Pog«Pd 4 i s Rayleigh number; M _ |a 0 Kfp measure the ratio of magnetic 
PK UK 1 o.p og (l + x) 

to gravitational forces, f = f ^Qd ^ j s ^ Baylor num ber, N = RM, = ^° — r- is magnetic thermal Rayleigh 
A I v J ' Mk(1 + x) 

number; M 3 = — -. c — measure the departure of linearity in the magnetic equation of state and values from one 

(i+x) 

(M 0 = %H Q ) higher values are possible for the usual equation of state. 
The dimensionless boundary conditions are 

w = 0, T = l,Dcp = 0 at z = Oandw = 0, T = 0,Dcp = 0at z = l. (19) 
Eliminating p from equation (16) we get 

f V2 " lV w + R ( ! + M ' ) V h° - RM^hDcp, + VtT § = 0 , < 20 > 

Pr dt J dz 



and 

I. Pv dtp 5z 

Where stands for z- component of the vorticity V H , is two-dimensional Laplacian operator on horizontal plane. 
Eliminating £, from equations (20) and (21), we get 

V 2 - — — 1 V 2 w + T A D 2 w + R(l + M, )V 2 H f V 2 - — — \ - PvM.V 2 f V 2 - — — Vm, = 0. 
Pvdt) V \ PrSt J 1 \ PrSt J 1 

III. NORMAL MODE ANALYSIS 

Analyzing the disturbances of normal modes and assume that the perturbation quantities are of the form 
[w, 0, <p ! ] = [W(z), 0(z),<D(z)]exp (ik x x + ik y y + nt) , 

where, k x? k y are wave numbers in x- and y- direction and n is growth rate of disturbances. 
Using equation (21), equations (20) and (17) - (18) becomes 

(("-••#-•>♦ 

W + (D 2 -a 2 -n)© = 0, 

D@-(d 2 -a 2 M 3 )3> = 0- (26) 
where i-* _ and a 2 = k 2 x + k 2 v is dimensionless the resultant wave number. 



jw-a 2 R(l + M^D 2 - 


-a'-^je + a'RM^D 2 -; 








PrJ 



(25) 



The boundary conditions of the problem in view of normal mode analysis are 

W = 0,D 2 W = 0, 0 = 0, DO = 0 at z = 0,l . (27) 

IV. METHOD OF SOLUTION 

The Galerkin weighted residuals method is used to obtain an approximate solute on to the system of equations (24) 
- (26) with the corresponding boundary conditions (27). In this method, the test functions are the same as the base (trial) 
functions. Accordingly W, ® and <D are taken as 

W = j^A p W p ,0 = j^B p 0 p ,DO = X C p D<I) p (28) 

Where A p , B p and C p are unknown coefficients, p =1, 2, 3...N and the base functions W p , 0 p and DO p are assumed in the 
following form for free-free boundaries are: 



W p = Cosp7i z , 0 p = Cosp 71 z, DOp = CospTt z , 



(29) 

Such that W p , 0 p and <D p satisfy the corresponding boundary conditions. Using expression for W, 0 and DO in equations 
(24) - (26) and multiplying first equation by W p second equation by 0 P and third by DO p and integrating in the limits from 
zero to unity, we obtain a set of 3N linear homogeneous equations in 3N unknown A p , B p and C p ; p =1,2, 3, ...N. For existing 
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of non trivial solution, the vanishing of the determinant of coefficients produces the characteristics equation of the system in 
term of Rayleigh number R. 

V. LINEAR STABILITY ANALYSIS 

We confined our analysis to the one term Galerkin approximation; for one term Galerkin approximation, we take 
N=l, we get the expression for Rayleigh number R as: 

(tT + a 2 M 3 Jn 2 +a 2 ] n 2 +a 2 +— +T A 7i 2 

R = ^ , FrJ s— ■ (30) 

a 2 (7i 2 + a 2 M 3 + a 2 M,M 3 J n 2 + a 2 + ^ j 

For neutral stability, the real part of n is zero. Hence we put n = ico, in equation (30), where ro is real and is dimensionless 
frequency, we get 



[n 2 + a 2 M 3 Jn 2 + a 2 )^7t 2 + i 


2 ico^ 


+ T A n 2 














a 2 {i^ +a 2 M 3 + a 2 M 1 M 3 ) 


[vl 2 +a : 











Equating real and imaginary parts of equation (31), we get 
a 2 R(n 2 +a 2 M 3 +a 2 M l M 3 \n 2 +a 2 



^M.+a^Mj^+a 2 ) (2 2 y 2 / J 2 2J 2 ^ 

(7T 2 +a 2 M 3 j v ' \Jr 2 Pr ) 

and 

r .'R(»' + .'M^.'M,M 1 ) _r ( , , y ^^.'y ^^ 0 

I Pr(7T 2 +a 2 M 3 ) ^ V ' Pr Pr 2 M 



From equation (36) it follows that co is real it is necessary that 
a 2 R(7i 2 +a 2 M 3 +a 2 M,M 3 Y7i 2 +a 2 ) / 2 2 Y ^ 2 

— ^ + a'Mj '*(«•+. ) + t a «. (34) 

Hence equation (32) gives the oscillatory stability provided that inequalities (34) hold, 
(a) Stationary Convection 

Consider the case of stationary convection i.e., co = 0, from equation (30), we have 

a [n +a M 3 +a M 1 M 3 ) 

From equation (35) we have 

5R _ 9R . 8R „ 

> 0, < 0 and < 0 . 

c7T A 5M 3 9M, 

Thus rotation has stabilizing effect while magnetization parameters have destabilizing effect on fluid layer. 
In the absence of rotation T A =0, the Rayleigh number R for steady onset is given by 

R _ (7T 2 + a 2 M 3 ^ 2 + a 2 ) 3 

a 2 (7i 2 +a 2 M 3 +a 2 M 1 M 3 ) 
This is the good agreement of the result as obtained by Finlayson (1970). 

In the absence of rotation T A =0 and magnetic parameters M!=M 3 =0, the Rayleigh number R for steady onset is given by 

K ^~?t 
a 

Consequently critical Rayleigh number is given by Rc = ^ K . 

4 

This is exactly the same the result as obtained by Chandrasekhar (1961) in the classical Benard problem. 

VI. CONCLUSIONS 

A linear stability analysis of thermal instability for ferrofluid in the presence of rotation is investigated for free- 
free boundary layer. Galerkin-type weighted residuals method is used for the stability analysis. The behavior of rotation 
and magnetization on the onset of convection analyzed. 
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The main conclusions are as follows: 

a. Rotation stabilizes the fluid layer while magnetization parameters destabilize the fluid in case of stationary 
convection. 

b. The 'principle of exchange of stabilities' is not valid for the problem. 

c. The oscillatory convection is possible if 

a 2 R(?i 2 +a 2 M, +a 2 M,Mj7T 2 +a 2 ) / 2 ,y _ 2 

i j~ -. — L r ^ ; -<[tz +a ) + T A n ■ 

7r 2 +a 2 Mj V ' 
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